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Abstract 

A hybrid grid system that combines the Chimera overset grid scheme and an unstruc- 
tured grid method is developed to study fluid flow and heat transfer problems. With 
the proposed method, the solid structural region, in which only the heat conduction is 
considered, can be easily represented using an unstructured grid method. As for the 
fluid flow region external to the solid material, the Chimera overset grid scheme has 
been shown to be very flexible and efficient in resolving complex configurations. The 
numerical analyses require the flow field solution and material thermal response to be 
obtained simultaneously. A continuous transfer of temperature and heat flux is speci- 
fied at the interface, which connects the solid structure and the fluid flow as an inte- 
gral system. Numerical results are compared with analytical and experimental data for 
a flat plate and a C3X cooled turbine cascade. A simplified drum-disk system is also 
simulated to show the effectiveness of this hybrid grid system. 



1. Introduction 


Prediction of a turbomachinery engine’s performance is a complex process entailing the 
iterative execution of aerodynamic, thermal, and structural analyses. Ultimately, engine perfor- 
mance is predicted from an aerodynamic analysis of the primary flowpath. However, the accuracy 
of the aerodynamic analysis is dependent on several factors, one of which is the definition of the 
flowpath geometry. Accurate determination of the flowpath geometry, which varies depending 
upon the operation conditions, requires the inclusion of aerodynamic, thermal and structural 
effects. Consequently, physical situations where heat transfer occurs between material and fluid 
flow with differing properties are commonly encountered in engineering practice, and often the 
geometries of interest are extremely irregular. 

In recent years, with advancements in performance of computing systems as well as 
improvements in numerical schemes, computational analyses of the Navier-Stokes equations for 
fluid flow and heat transfer has become increasingly feasible (Chima, 1985; Ni et al., 1986; 
Kwon, 1988; Boyle, 1991; Li et al., 1994). Although important progress has been made in solving 
the Navier-Stokes equations, much work is still needed to achieve robustness, accuracy, and effi- 
ciency. An issue which concerns the efficiency and accurancy of the design analysis for engine 
components is the grid system used for the numerical simulation. For instance, the turbine blades 
are often characterized by high turning geometry and large mean flow deviation from the axial 
direction, resulting in substantially skewed grids. The man power required to generate an ade- 
quate grid system for various engine components is usually expensive; the procedure is ineffi- 
cient. Another issue whcih also arises in the analysis is the thermal condition at the solid body 
surface. In conventional computations, the thermal conditions at the body surface are imposed 
with either a constant temperature condition or a prescribed heat flux condition when solving the 
Navier-Stokes equations. In reality, the temperature condition and the heat flux at the solid surface 
are both unknown a priori and must be determined as a part of the solution. As a result, a conju- 
gate analysis of the fluid flow and heat transfer must be undertaken to obtain physically realistic 
solutions. 
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Among various grid systems for numerical simulation, the Chimera overset grid method 
(Benek et al., 1985) has been shown to have the most flexibility and efficiency in dealing with 
very complex configurations (Buning et al., 1989). It consists of a set of overlapped structured 
grids which are independently generated and body-fitted, yielding a high quality grid readily 
accessible for efficient solution schemes. Therefore, the Chimera grid method can be utilized to 
facilitate the grid generation process, resulting in a more efficient design procedure. In the conju- 
gate problem, the fluid flow external to the solid body and the heat conduction inside the solid 
body must be solved in a coupled fashion. As is known in solid mechanics analysis, most struc- 
tural problem can be easily resolved by using an unstructured grid method to represent a rigid 
body. However, any pure-strain approach will not achieve the efficiency for fluid flow and heat 
transfer problems. Hence, a method that properly employs a hybrid of the Chimera and unstruc- 
tured grids may prove to be fruitful. A hybrid grid system using an unstructured grid to replace the 
arbitrarily overlapping region in the framework of the Chimera scheme, termed DRAGON grid, 
has been developed by Kao and Liou (1995). Their method conserves the numerical fluxes at the 
grid interface and shows some promising features in terms of the solution accuracy and the grid 
flexibility. In this paper we adopt their approach by using the unstructured grid to resolve the heat 
conduction inside the solid body while remaining with the Chimera overset grid scheme for the 
fluid flow region. 


2. Governing Equations and Numerical Algorithm 
The governing equations for the fluid flow are the 3-D, compressible Navier-Stokes equa- 
tions. The heat conduction equation is used to evaluate the temperature condition inside a solid 
material. 

2.1 Navier-Stokes Equations 

The time-dependent compressible Navier-Stokes equations, expressed in an integral form 
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over an arbitrary control volume Q are: 


J|?dv+ = 0 

n an 

where the conservative-variable vector U = (p, p u, pv, pw, p F) T . The inviscid flux is written as a 
sum of the convective and pressure fluxes: F = OV + P, where O = (p, p u, pv, pw, pH) T is 

A. V \ A 7* 

the vector containing convected variables and P = p (0, i,j, k, 0) . The specific total energy is 
E = e f . + |v| 2 /2 = H—p/p. We denote with an overhead arrow the vector quantities expressed 
in terms of Cartesian coordinates. 

Based on the cell-centered finite volume method, the governing equations are semi-dis- 
cretized. We use the new flux scheme AUSM + developed by Liou (1994) to express the numerical 
flux at the cell faces. The AUSM* scheme allows an exact capture of a normal shock by using a 
suitably chosen interface speed of sound, yields a smoother solution by way of including higher- 
order polynomials, and leads to faster convergence rate. 

The semi-discretized form, describing the rate of time change of U in £2 via the balance of 

A 

fluxes through all enclosing faces, S[, l = 1, .... LX, whether they be in the structured or unstruc- 
tured grid regions, can be cast as 



The last terms account for only normal components of the flux at the face, in terms of which the 
equations are written as: 

= 0 

where Fnl = Fl'*l and h l is the unit normal vector of 5/. The task is then to represent the 
numerical flux at the cell interface 5/, which straddles cells denoted by subscripts “L” and “R”. 
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The AUSM + scheme gives the numerical flux in the following expression. 

^ = Afl> + p m 

In this formula, a l is the speed of sound suitably defined at the interface to result in an exact reso- 
lution of a normal shock. Here M/ is the interface Mach number and P nl is the interface pressure 
in the surface normal direction. 

2.2 Heat Conduction Equation 

The conductive heat transfer in the solid structure is governed by the integral heat conduc- 
tion equation. 



Cl dCl 


where E = p CT denotes the potential energy of the solid material and q is the heat flux, 

q = -KVT 

Here is the thermal conductivity of the material, is the density and C is the specific heat. 

The above integral heat equation is then expressed using the finite volume method. The 
result is that the heat flux through a surface depends on the temperature in the neighboring cells 
near the surface. The change in temperature in a cell therefore can be obtained by integrating the 
heat fluxes throughout the surrounding cell faces. To be consistent with the flow solver, this sys- 
tem is solved iteratively, using a two-step predictor/corrector method on each time step. 

2.3 Time Integration 

The time integration scheme updates the conservative variables at the cell center. The 
present method originates from the Taylor series expansion in time, as was done by the Lax-Wen- 
droff scheme. Then a two-step scheme with second order time accuracy can be obtained. 
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Predictor. 


U* = lf + M-4-lf 1 

dt 

Corrector: 

U = U +A t-4-U 

dt 

U " +1 = ±(U n + U*) = ±{U +[U n + At-j i U^} 

It is noted that (1) the predictor step allows a full time step, and (2) like other 2-step integra- 
tion schemes, only two levels of storage are needed as if 1 is absorbed in forming the residual 
indicated in the squared bracket, (3) both predictor and corrector steps are identical, with no need 
of defining a midpoint for the corrector step, leading to simplification of coding and the complex- 
ity of evaluating the transport terms. 

The stability of the scheme is restricted by the CFL (Courant-Friedrichs-Lewy) number not 
exceeding unity. For steady state cases, local time stepping is used to accelerate conveigence. 

3. Grid Methods 

The Chimera composite grid method has been shown to be very flexible and efficient in 
generating body-fitted structured grids for complex configurations. In the solid mechanics com- 
munity, unstructured grid techniques are broadly used to represent structural domains. A combi- 
nation of these two grid systems is developed to tackle the conjugate heat transfer problems in this 
paper. 

3.1 Chimera Grid 

The Chimera scheme is a grid embedding technique which provides a conceptually simple 
method for domain decomposition. For instance, a major grid is generated about a main body ele- 
ment and minor grids are then overset on the major grid so as to resolve interesting features of the 
configuration. Usually the minor grids are overset on top of the major grid without requiring the 
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mesh boundaries to join in any apecial way. However, a common or overlap region is always 
required to provide the means of matching the solution across boundary interfaces. To increase 
the flexibility in the selection of subdomains, the Chimera scheme also allows a means of remov- 
ing regions of a mesh containing an embedded grid. That is, an embedded mesh introduces a 
“hole” into the mesh in which it is embedded. Typically, a hole is defined by a creation boundary 
which contains a surface or a group of surfaces. The purpose of a hole creation boundary is to 
identify points that are within this boundary. A mesh point is considered to be inside a hole cre- 
ation boundary if it is inside all surfaces that define the boundary. Figure 1 illustrates the connec- 
tions between composite overlapping grids, with hole points being blanked by a prescribed 
creation boundary. A practical application of the Chimera grid system for the complete Space 
Shuttle/Solid Rocket Booster geometry (Buning et al., 1989) is shown in Fig. 2. 

A flow solver must be modified to account for the use of multiple meshes and the “hole” in 
the grids. These hole points must be blanked or excluded from the flow field solution. The main 
change in the flow algorithm itself is the treatment of the hole boundaries. The blanked solutions 
are updated in the interpolation routine. With this approach, no special routine or logic tests are 
required to exclude the blanked points from the flow field solution. Details regarding the grid 
embedding technique can be found in Steger et al. (1987) and Suh et al. (1990). 

3.2 Unstructured Mesh 

The Delaunay triangulation scheme is applied to generate an unstructured grid in the solid 
material region. The steps for adopting the unstructured cells in the framework of the Chimera 
grid scheme are summarized below. 

1. Boundary nodes provided by the Chimera grid meshes are reordered according to their 
geometric coordinates (Barth, 1994). 

2. The Delaunay triangulation method is then performed to connect these boundary nodes 
based on the Bowyer algorithm (1981). 

3. In the unstructured grid, since their is no logical ordering of the cells and their neighbors, 
connectivity matrices containing cell-based as well as edge-based information are introduced. 
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Also, the present approach requires additional matrices to connect the structured and unstructured 
grids. The connectivity matrices used in the present 2D version are summarized as follows. 

(a) IEDGENODE(l:2, edges) - 2 nodes for each edge, 

(b) ICELLEDGE(1:3, cells) - 3 edges for each triangular cell, 

(c) IEDGECELL(1:2, edges) - 2 neighboring cells for each edge, 

(d) IEDGETYPE(edges) - edge type (type of boundary condition), 

(e) IEDGEFLUX(grids , edges) - edge number that connects structured and unstructured 
grids, where the parameter “grids” identifies which structured grid, 

(f) IFLUXINDX(grids, edges) - i-index of the structured cell that shares the edge fluxes 
with the unstructured cell, 

(g) IFLUXINDY(grids, edges) - ./-index of the structured cell that shares the edge fluxes 
with the unstructured cell. 

The first four types of connectivities are standard in unstructured grid codes, except that we 
assign “+” and values to IEDGETYPE to indicate whether the outward normal vector of a tri- 
angular cell face points into or out of the neighboring structured cell which shares the same face. 
This information is needed as the tempercture and heat flux conditions are enforced to be identical 
on this fase. 

Questions concerning triangular grid quality and improvements are not pursued for they are 
beyond the scope of the present paper. 

3.3 Data Communication Through Grid Interfaces 

In the Chimera method, communication between overlaid grids is achieved by interpolation 
of boundary values from the mesh or meshes in which the boundaries are contained. We have 
used the PEGSUS codes (Benek et al., 1986; Suhs et al., 1991) to make interconnections of sub- 
domains, define the hole regions, and supply points to facilitate communication among grids dur- 
ing the solution process. At present, the Chimera scheme employs bilinear interpolation that, as 
some simple experiments have shown (Mastin et al., 1984; Parthasarathy, 1984), is superior to 
Taylor series expansion. 
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In the current work, both the structured and unstructured solver are based on the cell cen- 
tered scheme in which quadrilateral and triangular cells are used repectively. Figure 3 shows the 
interfaces connecting both structured and unstructured grids. As noted, the heat flux, evaluated at 
the cell interface, is based on the conditions of neighboring cells. In the conjugate problem, the 
temperature and the heat flux values at the interface are continuously transferred according to the 
following relations: 


q w = -K :pT f = -xvr, 

Here the subscripts w,f and s denote the values at the solid body surface, the fluid boundary and 
the solid boundary, respectively. 


4. Test Results and Discussions 

Numerical results are compared with analytical and experimental data for a flat plate and a 
C3X cooled turbine cascade. A simplified drum-disk system is also simulated to show the effec- 
tiveness of the hybrid grid method. For the hybrid grids, both Chimera and unstructured grids can 
be displayed on the same plot using the Flow Analysis Software Toolkit (FAST; Walatka et al., 
1993) visualization package on an IRIS workstation. 

Case 1 - Conjugate Heat Transfer on a Flat Plate 

The first case is a flat plate with a cylindrical leading edge, as shown in Fig. 4. The upstream 
flow conditions are M = 03,Re l = 10 4 , Pr = 0.71 and T„ = 300 K. Here l = 1 m is the 
chord length of the plate and the wall thickness is 2.5 cm. Aluminum, with a density of 2700 kg/ 
m 3 , thermal conductivity of 211 W/m-°K and a specific heat of 900 W-s/kg-°K, is used as the 
plate material. The inner wall temperature is set to be 280°K throughout the computation. 

Under these conditions, the temperature along the outer wall is nearly constant and equal to 
the inner wall temperature. Therefore, the calculated Nusselt number can be compared with the 
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analytical values, which are given for a flat plate of zero thickness and constant wall temperature 
by Eckert et al. (1941). 


Nu = 0.33337 P?jRe~ x ■ l - 


The computed Nusselt number can be obtained by 


Nu 


<rj_ 

Kf(T w -T a ) 


where q is the heat flux, / is the plate length, T w is the wall temperature, T a is the adiabatic wall 
temperature, and xy is the thermal conductivity of the fluid. 

Figure 5 shows the hybrid grid system used for the flat plate case. The external flow field is 
first resolved using one structured mesh with grid dimensions of 168 x 50. The Delaunay triangu- 
lation method is then used to construct the unstructured grid to represent the flat plate. In this case, 
we use 530 nodes and 800 triangular cells for the unstructured grid region. The computed Nusselt 
numbers along the plate surface are compared with the analytical values and displayed in Fig. 6. 
The profiles are nearly identical everywhere except near the leading edge, where the analytical 
results are not valid. 

Case 2 - C3X Cooled Turbine Cascade 

In this case, the C3X cooled turbine cascade reported by Hylton et al. (1983) is simulated 
for the conjugate heat transfer solution. The blade geometry considered has been shown experi- 
mentally to exhibit a two dimensional response. The blade was instrumented with surface thermo- 
couples and pressures sensors. Data on the blade surface were acquired via a computer controlled 
data acquisition system, and reduced data are tabulated and reported in Hylton et al. (1983). 

A C3X test case (run 112) is selected and its set-up conditions are listed in Hylton et al. 
(1983). The coolant flow conditions and the experimentally determined external blade wall tem- 
perature, heat transfer coefficient, and pressure distributions are also described in that report. 

At the inlet, for run case 112, the Mach number is 0.17, the Reynolds number is 0.53xl0 6 , 
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the total temperature is 783 °K, and the total pressure is 6.22 KPa. The blade geometry possesses 
a high turning angle which causes difficulties in the grid generation process. A Chimera overset 
grid system, as displayed in Fig. 7, utilizes an O-type grid (size 139x51) wrapping around the tur- 
bine cascade, and an H-type grid (size 76x70) to cover the fluid flow path. This grid system can be 
applied with ease, not only to preserve high grid quality at the region near the blade body, but also 
to maintain the grid periodicity at the upper and lower boundaries. The Delaunay triangulation 
method is then used to generate the unstructured grid representing the inner region of the blade, 
with a grid size of 610 nodes and 1001 triangular cells. Note that the 10 cooling holes are also 
included to simulate the conjugate heat transfer solution with cooling effect. 

The Baldwin-Lomax algebraic turbulence model reported in Baldwin et al. (1978) is used to 
calculate the eddy viscosity. The effective viscosity is then given by 




tJ 


Here the subscripts / and t refer to the values in the laminar and turbulent flows. In the fluid flow 
and the solid body regions, the thermal conductivities remain unchanged throughout the computa- 
tions. The Prandtl numbers, Pr t and Pr t , are equal to 0.71 and 0.9 for laminar and turbulent 
flows, respectively. The solid material used for the blade body is stainless steel. 

Two types of blade boundary conditions are applied in this test case: (1) adiabatic wall, and 
(2) conjugate blade surface with internal cooling effect. Figures 8-9 illustrate the Mach number 
contours and the temperature distributions for the adiabatic and conjugate cases, respectively. The 
pressure distributions along the blade surface for both adiabatic wall and conjugate wall cases, as 
shown in Fig. 10, indicate close agreement with the experimental data. 

As the results plotted in Fig. 11 indicate the adiabatic wall temperature is significantly 
higher than the experimentally determined blade surface temperature. The prediction of the blade 
temperature is greatly improved when the conjugate heat transfer analysis with internal cooling 
effect is included. It is seen that the conjugate blade temperature displayed in Fig. 11 agrees well 
with the experimental data in both the pressure and the suction regions. However, the discrepancy 
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shown around the trailing edge of the blade indicates a higher-order turbulence model may be 
needed to improve the solution of the wake flow. 

Case 3 - Simplified Drum-Disk System 

A secondary flow system in the gas turbine engine is an example which includes aero, ther- 
mal and structural analyses. In this test case, a simplified drum-disk flow system is simulated 
using the present hybrid grid method. Rotating effects of the drum-disk system are not considered 
in the present calculation. 

As illustrated in Fig. 12, the Chimera overset grid scheme is used to represent the interior 
region between the drum and the disk. It is noted that the disk grid (88 x 21 ) is embedded on the 
drum grid (116x61). Consequently, those grid points inside a prescribed hole creation boundary 
have been blanked out from the drum grid. As also displayed in Fig. 12, an unstructured grid has 
been constructed to fill in the disk body, with a grid size of 195 nodes and 301 triangular cells. 

In this test case, two types of thermal boundary conditions are used at the disk-fluid conju- 
gate surface, namely the adiabatic wall and conjugate heat transfer wall conditions. Figures 13 
and 14 illustrate the pressure and Mach number contours, respectively. It is apparent from the 
plots that the effect of heat transfer has little influence on the fluid flow solution. 

The computed temperature distributions for the drum-disk system are laso displayed in Fig. 
15, indicating similarities in most contours in the fluid region except that some major differences 
are found near the disk surface. Comparisons of the temperature distributions along the disk body 
surface are ploted in Fig. 16, showing the average adiabatic wall temperature is much higher than 
the conjugate wall temperature. Subject to the adiabatic wall assumption, the heat flux is mainly 
convected by the fluid flow and its value vanishes at disk body surface, causing a large tempera- 
ture deviation. With the conjugate wall condition, the heat flux is conducted as well as convected 
within the fluid flow and solid body regions, promising a better thermal exchange at the disk sin- 
face. Note that the temperature variations along the surface of the disk body are rather small, and 
this trend is consistent with the design condition in the secondary flow system. 
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5. Conclusions 


In this paper we perform computations for conjugate heat transfer problems using a hybrid 
grid system. The fluid flow and solid structural regions are resolved by the Chimera overset grid 
scheme and unstructured grid method, respectively. In the computational process, the temperature 
and the heat flux at the cell interfaces between the fluid and solid regions are enforced to be con- 
tinuously transferred, with no artificial temperature distribution or prescribed heat transfer rate 
specified. The effects of heat conduction in the solid regions are shown to significantly affect the 
computed surface temperatures in comparison to those obtained from a traditional adiabatic wall 
assumption. With the proposed hybrid grid system, aerodynamic analysis for designing the engine 
components can now efficiently include aero, thermal and structural effects as one integral sys- 
tem. 
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Fig. 2 Chimera grid used for integrated Space Shuttle geometry. 
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Fluid Flow (Navier-Stokes) 



Solid Structure (Heat Conduction) 


Fig. 3 Heat flux at the cell face connecting the structured and unstructured grids. 
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Fig. 4 Conjugate heat transfer on a semi-infinite flat plate. 


17 



8 






Fig. 7 The Chimera/unstructured grid system used for C3X turbine cascade 
with 10 cooling holes. 


19 




J 


■ 0.0 0.2 0.4 0.6 0.8 1.0 

X/C 

Fig. 10 Pressure distributions along the blade surface. 



Fig. 11 Temperature distributions along the blade surface. 
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Fig. 13 Comparisons of the pressure contours for simplified drum-disk system. 

(a) adiabatic wall, (b) conjugate wall. (P max =1.44, Pmi n =0.34, AP=0.026) 
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(a) ad iabatic wall 



(b) conjugate wall 



Fig. 14 Comparisons of Mach number contours for the simplified drum-disk system, 
(a) adiabatic wall, (b) conjugate wall. (M max =0.7, M min =0, AM=0.03) 


(a) adiabatic wall 



(b) conjugate wall 



Fig. 15 Comparisons of the temperature contours for the simplified drum-disk system, 
(a) adiabatic wall, (b) conjugate wall. 

( T max =1 > T min=°- 6 ’ AT=0.01 - fluid region; 

T max=°- 7 ’ T min=°- 6 ’ AT=0.0004 - disk region) 
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